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Local compositions and structures of Zni_;tMg;tO alloys have been investigated by Raman and solid-state 
^^Zn/^^Mg nuclear magnetic resonance (NMR) spectroscopies, and by neutron pair-distribution-function (PDF) 
analyses. The E^^^ and E^^^^ Raman modes of Zni_;tMg;tO display Gaussian- and Lorentzian-type profiles, re- 
spectively At higher Mg substitutions, both modes become broader, while their peak positions shift in opposite 
directions. The evolution of Raman spectra from Zni_;tMg;tO solid solutions are discussed in terms of lattice de- 
formation associated with the distinct coordination preferences of Zn and Mg. Solid-state magic-angle-spinning 
(MAS) NMR studies suggest that the local electronic environments of ^^Zn in ZnO are only weakly modified by 
the 15 % substitution of Mg for Zn. ^^Mg MAS spectra of Zno.85Mgo.15O show an unusual upheld shift, demon- 
strating the prominent shielding ability of Zn in the nearby oxidic coordination sphere. Neutron PDF analyses 
of Zno.875Mgo.125O using a 2 X 2 X 1 supercell corresponding to Zn7Mg08 suggest that the mean local geome- 
try of Mg04 fragments concurs with previous density functional theory (DFT)-based structural relaxations of 
hexagonal wurtzite MgO. Mg04 tetrahedra are markedly compressed along their c-axes and are smaller in vol- 
ume than Zn04 units by ^ 6 %. Mg atoms in Zni_;tMg;tO have a shorter bond to the c- axial oxygen atom than 
to the three lateral oxygen atoms, which is distinct from the coordination of Zn. The precise structure, both 
local and average, of Zno.875Mgo.125O obtained from time-of-flight total neutron scattering supports the view 
that Mg- substitution in ZnO results in increased total spontaneous polarization. 

PACS numbers: 71.55.Gs, 77.22.Ej, 78.30.-j, 61.05.fm, 61.05.Qr 



I. INTRODUCTION 

Research in the area of polar semiconductor heterostruc- 
tures has been growing rapidly, driven in large part by inter- 
est in two-dimensional electron gas (2DEG) sy stems .^^'^'^1^ 
2DEGs are known to form at heteroj unction interfaces that 
bear polarization gradients. They can display extremely high 
electron mobilities, especially at low temperatures, owing to 
spatial confinement of carrier motions.^ Recent reports of 
2DEG behaviors in Gai_;,Al;,N/GaN and Zni_;,Mg;,0/ZnO 
heterostructures have great significance for the development 
of novel high-electron-mobility transistors (HEMTs)^^ and 
quantum Hall devices.^ 

2DEG structures are usually designed by interfacing a polar 
semiconductor with its less or more polar alloys in an epitax- 
ial manner. Since the quality of the 2DEG depends critically 
on interface perfection, as well as the polarization gradient 
at the heteroj unction, understanding compositional and struc- 
tural details of the parent and alloy semiconductors is an im- 
portant component in 2DEG design and fabrication. The evo- 
lution of atomic positions and cell parameters upon alloying 
can directly aff'ect the magnitude of the polarization gradient 
and the carrier density at the heteroj unction. 

Zni_;cMg;cO/ZnO is one of the more promising heterostruc- 
ture types for studies of 2DEGs, due to the large polariza- 
tion of ZnO, the relatively small lattice mismatch, and the 
large conduction band off'sets in the Zni-^^Mg^cO/^nO het- 



erointerface. Although 2DEG formation in Zni-^^Mg^cO/^nO 
heterostructures have been researched for some time, a clear 
understanding of the alloy structure of Zni_;cMg;cO is cur- 
rently lacking. Recently, we have studied composition- 
dependent changes in the crystal structures of poly crystalline 
Zni-x^gxO by synchrotron x-ray diffraction and Raman 
spectroscopy.l^El For the composition range < x < 0.15, 
we have shown that Mg-substitution modifies the aspect ratio 
of the hexagonal lattice through enhanced bond ionicity, and 
in parallel, decreases static polarization in the crystal, due to 
decreased internal distortion in the tetrahedral coordinations. 

Here, we conduct a detailed and more precise study of 
the local structure of Zni-x^g^O alloys using Raman and 
solid-state nuclear magnetic resonance (NMR), in conjunc- 
tion with neutron diff'raction techniques. Raman and NMR 
spectroscopy are useful probes for addressing the molecular 
compositions and structures of solid-solution systems. Peak 
shapes and widths of the Raman spectra reffect compositional 
ffuctuations and both short- and long-range order, whereas 
NMR is sensitive to the local environments around specific 
nuclei. Thus, both techniques can provide structural infor- 
mation on Zni_;cMg;cO that is complementary to diff'raction 
analyses. Here, we examine the details of and E^^^ Ra- 
man modes for Zni_;,Mg;,0 (x = 0, 0.05, 0.10, and 0.15). We 
have also used spin-echo magic- angle- spinning (MAS) ^^Zn 
and ^^Mg NMR to study separately zinc and/or magnesium 
species in ZnO and Zno.85Mgo.15O. These measurements have 
been correlated with average and local crystal structures of 
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Zni-JAgxO and ZnO, as established by Rietveld (diffraction 
space) and pair-distribution-function (PDF, real-space) analy- 
ses of time-of-flight neutron diffraction data. Compared with 
x-rays, neutron scattering provides much greater sensitivity to 
oxygen and Mg positions, as well as larger momentum trans- 
fer as measured by the larger maximum Q wave vector.^ By 
taking advantage of the increased data quality, we are able 
to isolate the geometry of Mg04 moieties stabilized in the 
wurtzite ZnO lattice. 



II. EXPERIMENTAL 

Poly crystalline powder samples of ZnO and Zni_;cMg;cO 
(x = 0.05, 0.10, 0.125, and 0.15) were prepared 
from oxalate precursors obtained by co-precipitation using 
Zn(CH3C02)2-2H20, Mg(N03)2-6H20, and H2C2O4, all of 
which had purities of 99.999 % from Aldrich. The two metal 
salts were dissolved together in deionized water and added to 
an oxalic acid solution in the ratio of [Zn^^]:[Mg^^]:[C204"] 
= (l-x):x:1.05. Upon mixing with oxalate, Zn^^ and Mg^^ 
immediately coprecipitated as white crystalline oxalate pow- 
ders, which were washed with deionized water and dried 
at 60°C for 4h and subsequently heated at 550°C for 20 h 
in air to decompose the oxalates to the oxides. Powder 
x-ray diffraction measurements confirmed the formation of 
Zni_;,Mg;,C204-2H20 after heating at 60°C, and Zni_;,Mg;,0 
after decomposing at 550°C. For the NMR experiments, sam- 
ples were also prepared at different conditions of temperature 
(900°C) and atmosphere (O2, N2). 

Raman spectra for ZnO and Zni_;cMg;cO (x = 0.05, 0.10, 
and 0.15) were acquired at room temperature using a Jobin 
Yvon-Horiba T64000 triple grating (1800gr/mm) spectrom- 
eter. The spectra of lightly compressed powders were 
recorded using micro-Raman sampling in air with 514.5nm 
and 675.5 nm excitations. The spectra reported here were av- 
erages of 10 acquisitions of 20 s integrations of the CCD de- 
tector with ^^lOmW of laser power focused onto the samples 
through a lOx objective. There was no evidence of degra- 
dation of the samples or associated changes in their spectra 
under these conditions. The spectrometer was calibrated us- 
ing a 520.7 cm"^ lattice mode of silicon. The Raman features 
of wurtzite E^^"^ and E^^^^ phonon modes were analyzed in de- 
tail to determine the peak shape, position, and width. Spectral 
background was removed following Shirle>^^, and the profile 
fitting was performed using pseudo-Voigi^ or Breit-Wigner- 
Fano functions. 

NMR measurements of the low-gyromagnetic-ratio nuclei 
^^Zn and ^^Mg were conducted at high (19.6Tesla) magnetic 
field strength at the National High Magnetic Field Labora- 
tory in Tallahassee, Florida. These investigations benefited 
from the enhanced sensitivity and improved resolution for 
the ^^Zn and ^^Mg nuclei that results from reduced second- 
order quadrupolar interactions, which scale inversely with the 
strength of the high applied magnetic field. The experiments 
were conducted at room temperature at 19.6 T (^H resonance 
frequency of 830 MHz), which for ^^Zn (/ = f , 4.1 % nat- 
ural abundance, ca. 1.5 receptivity relative to ^^C) and ^^Mg 



(/ = -|, 10.1 % natural abundance, ca. 0.7 receptivity relative 
to ^^C) corresponded to Larmor frequencies of 5 1 .88 MHz and 
50.76 MHz, respectively. All of the NMR spectra presented 
here were recorded on a single-resonance 4-mm probehead 
under MAS conditions at 10 kHz, using a Hahn-echo (i.e. |- 
T-Jt-T-acquisition) with the delays x set to one rotor period Xr. 
The ^^Zn and ^^Mg shifts were referenced to the bulk exter- 
nal secondary standards ZnSe (274 ppm relative to 1.0 M a£. 
Zn(N03)2) and MgO (26 ppm relative to 3.0 M aq. MgS04M 
respectively. Pulse lengths of 2 [is and 4 [is were used for the 
I and Jt pulses, respectively. A recycle delay of 1 s was used 
in each case, using 16,000 transients for the ^^Zn MAS spec- 
tra (ca. 5h each), and 160,000 transients for the ^^Mg MAS 
experiments (ca. 46 h each). Second-order quadrupolar MAS 
lineshapes were fitted using the program DMfit.^"^ 

Time-of-flight neutron diffraction data for ZnO and 
Zno.875Mgo.125O were obtained on the neutron powder diffrac- 
tometer NPDF at the Lujan Neutron Scattering Center at Los 
Alamos National Laboratory. For ease of structural modeling 
in the PDF analysis of Zni_;cMg;cO. we chose a composition of 
X = I and have used an appropriate 16-atom supercell model 
of ZnyMgOg. For each sample, 2 g of powder were packed 
in a vanadium can, and the data were collected for 3h at 
25°C using four detector banks located at 46°, 90°, 119°, and 
148°. Control runs, for the intensity corrections, employed an 
empty vanadium can, a vanadium piece, and air. The program 
PDFgetN was used to extract the PDF G(r) from the raw scat- 
tering data.CS First, the coherent scattering intensity I(Q) was 
obtained from raw data by the intensity corrections for con- 
tainer, background, and incident beam. Then, the I(Q) was 
converted to the structure factor 5(2) by the corrections for 
sample absorption, multiple scattering, and inelasticity. Fi- 
nally, the PDF G(r) was constructed by the Fourier transform 
of the reduced structure factor F(Q) = Q[S(Q) - 1]. In the 
Fourier transform of F(Q) to G(r), Q data were terminated at 
35 A" ^ PDF refinements were performed using the software 
program PDFgui.^^ The average crystal structures of ZnO and 
Zno.875Mgo.125O were determined by the Rietveld method us- 
ing the GSAS-EXPGUI software suite. 

III. RESULTS AND DISCUSSION 

A. Raman spectroscopy 

Room temperature Raman spectra for Zni-^^Mg^^O (x = 0, 
0.05, 0.10, and 0.15) are shown in Fig.J T] along with mode 
assignments for the observed peaks .1^^*^ The wurtzite lattice, 
with space group C^^ (Hermann-Mauguin symbol PS^mc), 
has four Raman-active phonon modes, Ai-\-Ei-\- 2E2. The two 
E2 modes are nonpolar, while the Ai and Ei modes are polar- 
ized along the z-axis and in the xj-plane, respectively.'EIlOlII] 
The polar modes are further split into longitudinal (LO) and 
transverse (TO) components due to the macroscopic elec- 
tric field associated with the LO modes. Raman spectra 
of the ZnO and Zni-x^gxO compounds are dominated by 
two intense peaks of (^ 98 cm'^) and E^'^^ (^ 436 cm'^) 
modes. Raman scattering efficiencies of individual modes in 
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FIG. 1 : Raman spectra recorded at room temperature for Zni__^Mg_^0 
(x = 0, 0.05, 0.10, and 0.15 from top to bottom). The excitation 
wavelength used was 514.5 nm. 
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FIG. 2: Line profile fittings for (a) E^^^ and (b) E^^^^ Raman peaks 
of Zni_;^Mg;»-0 (observed, gray circles; fit, solid lines; x = 0, 0.05, 
0.10, and 0.15 from top to bottom). Simplified lattice vibrations are 
depicted using large open circles for cations (Zn, Mg) and small gray 
circles for oxygen. 

ZnO are known to vary with the excitation energy.C^ With 
514.5 nm (2.41 eV) excitation, the highest Raman efficiencies 
are observed from E^^^^ and E^^"^ modes. However, the po- 
lar LO modes exhibit a strong resonance effect as the excita- 
tion energy approaches the electronic transition energies. In 
cases when ultraviolet lasers are used for excitation, the Ra- 
man spectra of ZnO or Zni_;cMg;cO are dominated by the sig- 
nals from LO modes. 

In order to quantify composition-dependent changes in the 



coq (cm ^) r (cm ^) 
437.0 I 1 I 




I \ I \ I 5 I I I I \ 

0.00 0.05 0.10 0.15 ' 0.00 0.05 0.10 0.15 



X in Zn^_^Mg^O 

FIG. 3: (a) Raman shifts coq and (b) linewidths F for the E^"^ (open 
circles) and E^^^^ (filled circles) modes of Zni_;(;Mg;(;0. 



Raman spectra of Zni_;cMg;cO. the E^^^ and ^2^^^ line pro- 
files were analyzed by least- squares fitting with standard peak 
functions. From the fits, peak position (coq) and the full- 
width- at-half-maximum (FWHM, F) linewidths were deter- 
mined. The profile fittings of both E2 modes are shown in 
Fig.[2j and the evolution of peak parameters coq and F with 
Mg-substitution x are plotted in Fig. [3] Consistent with the 
previous study,^ the E^^^^ peaks of Zni_;cMg;cO are best rep- 
resented by Lorentzian Breit-Wigner-Fano (BWF) lineshapes 
with some asymmetry. However, the E^^"^ line profiles do not 
fit as well to the Lorentzian-type functions. Instead, the E^^"^ 
lines of Zni_;cMg;cO are found to be pseudo-Voigt-type with 
predominantly Gaussian components of 97%, 83%, 77%, and 
81 % for X = 0, 0.05, 0.10, and 0.15, respectively. 

As shown in Fig. [3^, the two E2 modes of Zni_;cMg;cO ex- 
hibit distinct dependences of phonon energy on the compo- 
sition. With increased concentration of magnesium, the E^^"^ 

mode shows a blueshift, while the E^^^^ mode exhibits a red- 
shift. For explaining these opposing trends, the vibrational 
eigenvectors of the wurtzite E2 modes need to be considered. 
As shown in Fig.[2| both E^^"^ and E^^^^ modes are associated 
with atomic motions in the aZ^-plane. The lower energy branch 
corresponds mainly (^^85 %) to the vibrations of heavier com- 
ponents (cations, in case of Zni_;cMg;cO), and conversely the 
higher energy one corresponds mainly to those of lighter com- 
ponents (oxygen). Consequently, the E^^"^ mode energy of 
Zni_;fMg;cO is explicitly aflTected by the cationic substitution, 
according to the reduced mass eflTect. Replacement of Zn with 
Mg will decrease the reduced mass of the oscillator and in 
turn increase the phonon energy. By comparison, the change 
of cation mass should have less influence on the E^^^^ mode 
energy, and its redshift is adequately attributed to the phonon 
softening caused by the in-plane lattice expansion. As previ- 
ously reported, the lattice constant a of Zni-x^g^O increases 
monotonically with x,^ which appears to account for the ob- 
served E^^^^ mode behavior. 



Peak shapes and linewidths of Raman spectra are dictated 
by the underlying Hne-broadening mechanisms. Gaussian 
hne broadening is intrinsic to Raman spectra and originates 
from instrumental resolution. Also, if the sample being ex- 
amined contains inhomogeneous components with different 
phonon frequencies, the resulting Raman lineshape would 
have a Gaussian character, reflecting the statistical nature of 
the spectrum. On the other hand, Lorentzian line broaden- 
ing is mostly attributed to finite phonon lifetimes (0, accord- 
ing to the energy-time uncertainty relationship Ylfi = l/t. 
The phonon lifetime shortening, which will cause linewidth 
broadening, can occur via two independent phonon decay 
mechanisms; by anharmonic decay into other Brillouin zone 
phonons or by phonon disruption at crystal defects. The lat- 
ter, involving lattice imperfections, are quite common for al- 
loy and solid- solution systems. Specifically for Zni-^^MgjcO, 
it is expected that Mg- substitution will substantially reduce 
the size of ordered domains, and thereby block the long-range 
propagation of zone center phonons. 

As noted above, the E^^^ and ^2^^^ lines of ZnO are respec- 
tively fitted with Gaussian and Lorentzian functions. It in- 
dicates that the E^^^^ mode of ZnO has anharmonic phonon 
decay channels, whereas the E^^"^ phonon does not. The 

Lorentzian line broadening behavior of E^^^^ mode, in Fig.jsJ), 
further manifests that the Mg-substitution increasingly popu- 
lates the phonon field of Zni_;cMg;cO with defect centers. The 
£"2°^ mode also undergoes a gradual line broadening with x, 
and its peak shape changes upon Mg-substitution. While un- 
substituted ZnO displays a nearly pure Gaussian peak for the 
£^2°^ mode, the peaks from Zni_;cMg;cO (x > 0) samples are 
found to contain ^20% of Lorentzian components. Therefore 
we conclude that, for the Zni-^^Mg^^O system, the same line 
broadening mechanism applies to both E^^^ and E^^^^ modes: 
phonon lifetime shortening by increased crystal defects upon 
Mg substitution. 



B. Solid-state ^^Zn and ^^Mg NMR spectroscopy 

A diff'erent view of the disorder induced by Mg-substitution 
in ZnO materials can be obtained from solid-state NMR, us- 
ing the NMR-active isotopes ^^Mg (spin / = -|) and ^^Zn 
(/ = |). These nuclei are challenging to observe, because their 
low natural abundances and low gyromagnetic ratios (see Ex- 
perimental) yield low signal sensitivities, and because second- 
order quadrupolar interactions associated with nuclear spins 
|/| > I can lead to a dramatic broadening of their NMR spec- 
tra. Second-order quadrupolar eff'ects, however, scale with the 
inverse of the static magnetic field ^o, and are therefore miti- 
gated at high magnetic fields, which also increase signal sen- 
sitivity. 

SoHd-state ^^Zn MAS NMR spectra acquired at 19.6 T of 
ZnO and Zno.85Mgo.15O are shown in Fig.|4] The materials 
were prepared by heating at 900° C in N2 and yield similar 
spectra that clearly exhibit scaled second-order quadrupolar 
powder MAS lineshape features. Associated values for the 
isotropic ^^Zn chemical shift, 5iso, the quadrupolar coupling 
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FIG. 4: Solid-state spin-echo ^^Zn MAS NMR spectra acquired 
at 19.6 T, room temperature, and 10 kHz MAS: (a) ZnO and (b) 
Zno.85Mgo.15O powders, both prepared at 900°C in N2. Accompany- 
ing the ZnO spectrum in (a) is a fit (dashed line) of the experimental 
lineshape to a simulated second-order quadrupolar powder pattern. 



constant, vq, and the asymmetry parameter, t]q, can be de- 
termined by fitting each resolved second-order quadrupolar 
MAS lineshape. A simulated powder pattern is shown as 
the dashed line in Fig.|4^ accompanying the ^^Zn MAS spec- 
trum of ZnO and corresponding to the calculated values 5iso = 
236 ppm, Vq = 378 kHz, and t]q = 0.1. Several discrepancies 
between the calculated and experimental powder patterns may 
be due to modest anisotropy of the powder particles and/or to 
uneven excitation of the diff'erent crystallite orientations by 
the radiofrequency pulses. The powder pattern and fit param- 
eters are nevertheless consistent with a single crystallographi- 
cally distinct and relatively symmetric ^''^Zn site in the wurtzite 
ZnO structure. The ^^Zn MAS spectrum of Zno.85Mgo.15O in 
Fig. 4b exhibits a similar lineshape and features, which are 
broadened by the presumably random incorporation of Mg 
atoms within the Zn-rich wurtzite lattice. 

The close similarity between the ^^Zn MAS NMR spectra 
of ZnO (Fig.|4^) and Zno.85Mgo.15O (Fig.|4]3) indicates that the 
substitution of 15 % of Zn atoms by Mg has little influence on 
the local electronic environments of a majority of the ^^Zn 
nuclei. Assuming that Mg atoms are distributed randomly, the 
probability of having no Mg atom in the second coordination 
sphere of a given ^^Zn nucleus is only 14 % (the first coordi- 
nation sphere consisting of four O atoms). Significant eff'ects 
on the ^^Zn spectrum would be expected if the quadrupolar in- 
teractions were substantially perturbed by the substitution of 
Mg atoms for 1, 2, 3, 4 (probabilities of 30 %, 29 %, 17 %, and 
7 %, respectively) or more of the 12 Zn atoms in the second 
coordination sphere. However, despite the high probability of 
modifications of the local chemical environment of the ^^Zn 
nuclei upon 15 % Mg-substitution, dramatic changes of the 
^^Zn spectrum are not observed, rather only modest broaden- 
ing that reflects increased disorder. Under these conditions, 
the difl'erences induced in ZnO or Zno.85Mgo.15O by changing 
the synthesis atmosphere (e.g., air, O2, or N2) and temperature 
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FIG. 5: Spin-echo ^^Mg MAS NMR spectra, recorded at 19.6 T, of 
(a) poly crystalline MgO powder (external reference for ^^Mg NMR 
shift), (b) Zno.85Mgo.15O powder prepared at 550° C in air, and (c) 
Zno.85Mgo.15O powder prepared at 900° C in N2. 



(over the range 550°C to 900° C) were too subtle to be clearly 
reflected in the corresponding ^^Zn MAS NMR spectra (data 
not shown), which are similar to those in Fig. [4^, b. Thus, it 
appears that ^^Zn NMR measurements of these ZnO-related 
materials are sensitive principally to substitutions that result 
in modifications of the first coordination sphere of the zinc 
nuclei (substitution of the anionic species), or to substitutions 
resulting in larger lattice distortions or increased electronic 
disorder than observed here.^^ 

Solid-state ^^Mg MAS NMR measurements complement 
the ^^Zn NMR results by being sensitive to the structural and 
electronic perturbations experienced directly by Mg atom in- 
corporated into ZnO lattices. The ^^Mg MAS NMR measure- 
ments were similarly conducted at 19.6 T and the resulting 
spectra are shown in Fig. [5] for poly crystalline MgO and dif- 
ferent Zno.85Mgo.15O powders prepared at 550°C in air and 
at 900°C in N2. The ^^Mg MAS spectrum in Fig.js^ of 
MgO shows a narrow well-defined ^^Mg peak that reflects 
the highly symmetric coordination environment of the single 
type of Mg site in its rock-salt structure. By comparison, 
the ^^Mg MAS spectra of the difl'erent Zno.85Mgo.15O powders 
contain broad and relatively unstructured lineshapes centered 
at ca. -30 ppm to -40 ppm that reflect broad distributions 
of signal intensity and thus broad distributions of local ^^Mg 
environments. The first Mg coordination sphere (Mg04) in 
Zno.85Mgo.15O is expected to be significantly distorted from 
its regular tetrahedral geometry, and the second coordination 
sphere, composed of Mg(Zn,Mg)i2, is expected to have an 
even larger number of diff'erent local configurations, bond dis- 
tances, and/or bond angles. The ^^Mg MAS NMR spectra 
appear to be more sensitive than the ^^Zn NMR results to 
local material environments in Zno.85Mgo.15O. However, the 
low signal-to-noise and broad lines observed for natural abun- 



dance ^^Mg in Zno.85Mgo.15O, even at 19.6 T, preclude a de- 
tailed and reliable analysis of the isotropic shift distributions 
and second-order quadrupolar broadening that would be nec- 
essary to extract directly quantitative information on nature of 
the Mg site distributions and/or their disorder. It is notewor- 
thy that the Zno.85Mgo.15O powders do not display any spec- 
tral feature corresponding to bulk MgO, thereby ensuring the 
phase purity of Zni-^^MgjfO (x < 0.15) solid solutions, within 
the sensitivity limit of the ^^Mg NMR measurements.*^ 

In common for ^^Si, ^^Al, and ^^Mg, an increase of the co- 
ordination number for a given cationic center counts for en- 
hanced local shielding, resulting in an upfield displacement of 
the signal(s) to lower shift values. ^^^^ The spinel MgAl204 
containing Mg04 fragments exhibits a ^^Mg shift of 52 ppm, 
compared to 26 ppm for MgO which is comprised of Mg06 lo- 
cal units. By comparison, clay minerals in which Mg atoms 
are coordinated by six (0,0H) ligands yield ^^Mg shifts in the 
range of to -lOOppm.^^ The wurtzite oxide Zno.85Mgo.15O 
has a 4-coordinated geometry of Mg, similar to MgAl204, but 
the former clearly shows a more upfield- shifted ^^Mg reso- 
nance at -30 40 ppm. We infer that ^^Mg MAS NMR 

signals are considerably influenced by the ligand identity be- 
yond the first coordination shell, and that Zn has a far stronger 
shielding contribution than Al, Si, or alkaline-earths, thus ac- 
counting for the observed resonance values. 

C. Neutron diffraction 

In previous studies, the crystal structure of ZnO has 
been reported [space group PS^mc, Zn at (||0) and O at 
(||w)] with lattice constants of a = 3.2427-3.2501 A and 
c = 5 . 1 948 - 5 .207 1 A, and atomic position parameter u = 
0.381-0.3 826! ^^ ' ^^ ' ^^ ' ^ ^ ' ^^ ' ^ ^ In our previous synchrotron x-ray 
study on poly crystalline ZnO, we obtained <2 = 3.2503 A, c = 
5.2072 A, and u = 0.3829.^ Here, we reexamine the crystal 
structures of poly crystalline ZnO and Zno.875Mgo.125O utiliz- 
ing neutron scattering, which has several advantages over x- 
rays. The coherent neutron scattering length of O (5.804 fm) 
is comparable to those of Zn (5.680 fm) and Mg (5.375 fm),^ 
as distinct from x-ray diffraction, where scattering is strongly 
weighted by Zn. Moreover, neutron scattering can provide 
high Q diffraction data with much less attenuation than x-rays, 
due to there being no fall-off' with Q in the form factor. 

Rietveld refinements of ZnO and Zno.875Mgo.125O were car- 
ried out on four histograms of data collected at detector loca- 
tions of 46°, 90°, 119°, and 148°. The highest-Q data from the 
148° detector covers the range of J > 0.25 A (Q < 25 A"^). 
The model used for the refinement in space group PS^mc has 
Zn/Mg at the 2b Wyckoff' position (||0) and O also at 2b posi- 
tion (||w). For Zno.875Mgo.125O, Zn and Mg were statistically 
distributed over the common site with fractional occupancies 
fixed to the respective compositions. The refinement con- 
verged with refiability factors of 7?wp = 2.57 %,7?p = 1.80%, 
and 



2.71 for ZnO, and R. 



wp 



2.89%,7?n 



2.04 %, 



and X = 3.56 for Zno.875Mgo.125O. The Rietveld refinement 
profiles and the structural parameters are given in Fig. [6] and 
Table[I| respectively. The present refinement results agree 
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FIG. 6: Rietveld refinement of time-of-flight powder neutron diffrac- 
tion profiles of (a) ZnO and (b) Zno.875Mgo.125O. Calculated patterns 
(gray lines) are superimposed on observed data (open circles), with 
Bragg positions and diff'erence profiles at the bottom. 
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FIG. 7: Average wurtzite parameters c/a and u of ZnO and 
Zni_;<^Mg;<^0 determined by Rietveld refinements of neutron (dia- 
monds) and synchrotron x-ray (squares, ref. ^) diffraction patterns. 
Note the significantly smaller error bars that emerge from the neu- 
tron study. The inset shows the data observed from existing wurtzites 
(see ref.— for details). 
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FIG. 8: Comparison of high-g scattering data from ZnO measured 
using synchrotron x-ray and pulsed-neutron radiation: (a,b) F(Q) = 
QIS (2) - 1], and (c,d) G(r). In the Fourier transform of F(Q) to G(r), 
Q data were terminated at 28 for x-rays, and 35 for neutrons. 



well with the findings from our previous synchrotron x-ray 
study.ElBoth studies indicate that Mg-substitution expands the 
a-parameter, compresses the c-parameter, and decreases the 
oxygen position parameter u, as summed up in Fig. [7] The 
two data sets from x-ray and neutron reveal an identical trend 
in the structural evolution in terms of a, c, and u. 



TABLE I: Structural parameters for ZnO and Zno.875Mgo.125O deter- 
mined from the Rietveld refinement of time-of-flight neutron diffrac- 
tion pattern. 





ZnO 


Zno.875Mgo.125O 


V(A3) 


47.603(1) 


47.578(1) 


a (A) 


3.24945(1) 


3.25058(2) 


c(A) 


5.20574(3) 


5.19771(3) 


c/a 


1.60204 


1.59940 


u 


0.38275(5) 


0.38214(6) 


t/i,o(Zn/Mg) (A2) 


0.00648(8) 


0.0065(2) 


t/iso(0) (A2) 


0.0101(1) 


0.0087(2) 


dM-o (A) 


(Ix) 1.9925(3) 


(Ix) 1.9868(3) 




(3x) 1.9729(1) 


(3x) 1.9743(1) 


^O-M-O (°) 


108.022(7) 


108.081(9) 




110.881(7) 


110.825(9) 




FIG. 9: (a) Probability for a cation in Zni_^Mg^O to have 0, 1,2, 
3, 4 or 5 Mg neighbors in the 12-membered coordination shell (r ^ 
3.25 A), and (b) mean separation between the nearest Mg-Mg pairs. 



For the composition range up to x = 0.15, the Zni_;cMg;cO 
solid solutions remain isostructural with wurtzite ZnO, with 
both Zn and Mg tetrahedrally coordinated by O atoms. How- 
ever, since the two cations Zn and Mg have clearly dis- 
tinct crystal chemistries in oxidic environments, it is con- 
jectured that they may have distinct local geometries within 
Zni_;cMg;cO lattices. In order to directly address this question, 
we have carried out neutron PDF analyses of Zno.875Mgo.125O 
and ZnO. FigurejS] shows the PDF data for ZnO, comparing 
the present neutron data with the x-ray data previously mea- 
sured at Beamline 11-ID-B of the Advanced Photon Source at 
Argonne National Laboratory. As can be seen from Figs. [8^ 
and[8]3, the neutron study provides scattering information over 
a significantly higher Q range and also with higher resolution. 
Correspondingly, the neutron PDF G(r) can be obtained for 
far wider r-ranges than when using x-rays (Fig.jSj:). Given 
this, neutrons are expected to enable better description of the 
Zni-xMgxO alloy structure(s). 

In order to obtain an appropriate real-space model for the 
PDF analysis of Zno.875Mgo.125O, we estimate as follows, 
the impact of Mg- substitution on the chemical environment 
within the ZnO lattice. Cations in the wurtzite lattice form 
a hexagonal close-packed sublattice, where each cation has 
12 nearest neighbors. At the Mg- substitution level of x = 
0.125, it is estimated that ^^80% of the cations have at least 
one Mg neighbor in the cation sublattice (Fig. [9]). Each cation 
(whether Zn or Mg) has as neighbors, 1.5 Mg atoms on av- 
erage, and the mean Mg-Mg distance is as short as 6.5 A. 
Therefore the coordination geometry of Mg04 will serve a 
substantial structural factor in the PDF of Zno.875Mgo.125O and 
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FIG. 10: (Color online) ZnvMgOg supercell model structure used for 
neutron PDF analyses of Zno.875Mgo.125O, and coordination environ- 
ments for three distinct Zn sites therein. In the top part, unlabeled 
oxygen atoms belong to the 03 type (see text). 



needs to be treated with independent parameters in the PDF 
refinement. 

On the basis of the above considerations, the structural 
model for Zno.875Mgo.125O was designed as follows. To dis- 
tinguish the Mg and Zn atoms in the unit cell, while main- 
taining the composition, a 2x2x1 supercell of the primi- 
tive wurtzite structure was selected, with the cell parameters 
a = b ^ e.sKc ^ 5.205 A, a = |3 = 90°, y = 120°, and 
the unit content Zn7Mg08. In the average structure scheme, 
in which the atomic positions are set as (||0) for cations and 
(||w) for anions, the external and internal geometries of all 
the tetrahedra are uniformly defined by the c/a ratio and u, 
respectively. However in the supercell model, we allowed the 
Mg04 unit to have its own coordination geometry, irrespec- 
tive of the Zn04 geometry. As a result, the eight O atoms in 
the supercell unit are divided into three groups; one O atom 
bonded to Mg along the c-axis (01), three O atoms bonded to 
Mg laterally (02), and the remaining four O atoms (03). The 
seven Zn atoms are also grouped according to their proximi- 
ties to the Mg atom; three Zn atoms having an apical 02 oxy- 
gen atom (Znl), three Zn atoms bonded to 02 atoms laterally 
(Zn2), and one Zn bonded to only 03 atoms (Zn3). Figure[TO| 
illustrates the supercell configuration and the three diff'erent 
types of Zn04 geometries. 

The PDF of ZnO was analyzed using a simple wurtzite 
cell (Zn202) obeying the symmetry requirements of the 
PS^mc space group. However, for the supercell refinement 
of Zno.875Mgo.125O, several symmetry constraints were lifted 
from the average structure description. The z-coordinates of 
the 01, 02, and 03 types of oxygen atoms were indepen- 
dently varied, and the z-coordinate of Mg was also refined. 
We also attempted in-plane displacements of the 02 atoms. 
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FIG. 11: Neutron PDF refinements for (a) ZnO and (b) 
Zno.875Mgo.125O, parts of the 1.2-22 A fittings. Calculated PDFs 
(gray lines) are superimposed on experimental data (open circles), 
with the diff'erences shown at the bottom of each. 



but found that these were unstable in the refinement. There- 
fore, the refinement for Zno.875Mgo.125O included three more 
position parameters than those for ZnO. For both ZnO and 
Zno.875Mgo.125O, the refinements of G(r) used the variables of 
lattice constants, atomic position parameters, isotropic tem- 
perature factors, structure scale factor, and peak sharpening 
coefi&cients. 



Figure 1 1 shows the neutron PDF refinements for ZnO and 
Zno.875Mgo.125O, carried out in the r range of 1.2 A to 22 A. 
The refinement parameters are summarized in Table|ll| For 
both structures, the refinements were achieved with satisfac- 
torily low /?w's of 7.9 % (ZnO) and 8.3 % (Zno.875Mgo.125O). 
In comparison with the Rietveld results, the PDF-refined lat- 
tice constants are systematically larger by ^^0.06%, but good 
agreements are observed for the cla ratios. Also, the oxy- 
gen atom position in ZnO is reasonably well reproduced from 
the PDF and Rietveld refinements. An interesting finding 
emerges from the atomic coordinates of Mg and its adjacent 
oxygen atoms (01 and 02) in Zno.875Mgo.125O. As found in 
Tablejnj the three groups of O atoms have clearly distinct z- 
coordinates from one another. Both the 01 and 02 groups are 
vertically shifted towards their respective Mg atoms. Mg it- 
self also moves towards the basal plane formed of 02 atoms. 
The resulting Mg04 unit has a squashed tetrahedral geome- 
try (point group Csv), with a much shorter height than those 
of Zn04 tetrahedra. However the 03 atoms, which are not di- 
rectly bonded to Mg, have similar z-coordinates to those found 



2.04 



^ 2.00 

CD 
O 

."S 1.96 



1.92 




Mg Znl Zn2 Zn3 




3 

> 

00 



FIG. 12: (a) Mg-0 and Zn-0 bond distances in Zno.875Mgo.125O, 
as analyzed by the neutron PDF refinement of a ZnyMgOg supercell 
model (c- axial bonds, circles; lateral bonds, diamonds), and (b) bond 
valence sums for each cation site. Horizontal lines indicate the Zn-0 
bond distances (in the upper panel) and BVS for Zn (in the lower 
panel) in ZnO. 



in ZnO. Therefore, the (Zn3)(03)4 tetrahedra retain the unper- 
turbed geometry of ZnO. The other types of Zn04 tetrahedra 
share corners with Mg04 and are likely to undergo intermedi- 
ate distortion. 

To examine the presence of short-range structural character- 
istics, the PDF refinements were carried out from 1.2 A to dif- 
ferent values of rmax between 10-30 A. The lattice constants 
and position parameters showed insignificant variations with 



and the values had shallow minima at 



22 A for 



both ZnO and Zno.875Mgo.125O. The Mg-0 and Zn-0 bond 



distances plotted in Fig. 12 represent the averaged results of 
multiple refinements over 1.2-rniax (A), where rmax was in- 
creased from 10 to 30 A in 2 A steps. From Fig. 12 1, we 



TABLE II: Structural parameters from the neutron PDF refinements 
of ZnO and Zno.875Mgo.125O, over the r range of 1.2-22 A. 



ZnO 



Zno.875Mgo.125O 



a (A) 
c(A) 
cla 
z{0) 



t/i,o(Zn/Mg) (A2) 
t/iso(0) (A2) 



3.2514(3) 3.2527(4) X 2 

5.2089(6) 5.2015(8) 

1.6020 1.5991 

0.3824(4)/0.8824(4) 0.3825(8)/0.8825(8) (03) 
0.386(3) (02) 
0.858(5) (01) 
0.486(8) 

0.0086(8) 0.0063(6) 

0.0070(6) 0.0081(8) 
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find a meaningful distinction between the Mg04 and Zn04 
bonding structures. The Mg-0 bonds are shorter than the 
Zn-0 bonds, and the former set consists of one short and 
three longer bonds, in contrast to the latter. These results 
are consistent with the greater sensitivity of the ^^Mg NMR 
spectra (Fig.jSj), c) to Mg substitution into ZnO-rich lattices, 
compared to the ^^Zn NMR results (Fig.|4]). For Mg and the 
three Zn atom types in Zno.875Mgo.125O, bond valence sums 
(Bvs:P are evaluated, as shown in Fig.[T2|). The valence 
sums for all of the cation and anion sites in Zno.875Mgo.125O 
are within a reasonable range around the ideal magnitude of 
2 valence units (v.u.). The relatively underbonded situation 
for Znl atoms and the slight overbonding of the Mg atom are 
the results of a rather drastic shift of the bridging 01 atom. 
The global instability index (GII)^^ of Zno.875Mgo.125O is cal- 
culated as 0.116v.u., which is significantly larger than that 
of ZnO (0.067 v.u.). The GII value is often used as a mea- 
sure of the residual bond strain, and is known not to exc eed 
0.2 v.u. for ordinary structures in the standard state.E^^It is 
worth mentioning, in light of the above GII considerations, 
that the increased microstrain upon Mg- substitution results in 
characteristic optical and x-ray line broadening behavior of 
Zni_;cMgjcO that we have previously reported.'^^*^ 

The solid solution of Zni-x^g^O has the thermodynamic 
solubility limit of x ^ 15 % on the ZnO-rich side. Obviously 
the distinct coordination preferences of Mg and Zn prevent 
the formation of a continuous solid solution across the entire 
composition range, x. While Mg has a point symmetry of Oh 
in its binary oxide MgO, wurtzite ZnO possesses only a Csv 
environment. Here, we rationalize the neutron PDF result for 
the Mg04 tetrahedral geometry, in several ways. Both bond 
valence^^ and ionic radii"^^ compilations are useful to explain 
the smaller tetrahedral volume of Mg04, as compared with 
Zn04. The bond valence parameter of Mg-0 (Rq = 0.693 A) 
is smaller than that of Zn-0 (Rq = 0.704 A), implying that 
in general Mg-0 bonds are shorter than Zn-0 bonds, for 
the same coordination numbers. Similarly, the 4-coordinate 
ionic radius of Mg^^ (0.57 A) is smaller than that of Zn^^ 
(0.60 A). In order for the Mg04 tetrahedron to have a smaller 
volume than Zn04, either the tetrahedral height or the base 
area should shrink, or an isotropic volume change could occur. 
From the PDF analysis, the c parameter around Mg showed a 
clear decrease, whereas the a parameter appears to be nearly 
unchanged. The c-axial parameter contraction of Mg04 unit 
can be deduced also from the composition-dependent changes 
of the average structure of Zni_;cMg;cO. 

The variable u is the only position parameter in the wurtzite 
structure. It corresponds to the ratio of the apical bond dis- 
tance to the c length, Japicai/<^. or <iapicai/2/i, where h is the 
height of the tetrahedron. For Zno.875Mgo.125O, u(Mg) is de- 
termined to be 0.395, a value that finds no correspondence 
to the experimental structures of pure wurtzites. Existing 
wurtzites have u values within the range of 0.374 to 0.383. 
Such a high u found for Mg04 is closely linked with its ex- 
tremely small aspect ratio (2h/a = 1.509, as compared to the 
ideal value 1.633). To compare the potential fields within the 
tetrahedral spaces of Mg04 and Zn04, the valence sums of 




0.50 



FIG. 13: w-dependent bond valence sums for Mg and Zn in the tetra- 
hedral spaces of their oxides. The dimensions of Mg04 and Zn04 
tetrahedra are from the neutron PDF refinements of Zno.875Mgo.125O 
and ZnO, respectively. For u near 0.5, a fifth oxygen atom was added 
to the coordination shell to account for the pseudo-bipyramidal con- 
figuration. 



ZnO, the tetrahedral cavity is slightly oversized for Zn, and the 
BVS(Zn) has a minimum of 1.9 v.u. at w ^ 0.382. The poten- 
tial minimum suggested by BVS(Zn) is in a good agreement 
with the experimental u for ZnO. However, the Mg04 shell in 
Zni_;cMg;cO is rather small for Mg, resulting in a global over- 
bonding situation. The BVS(Mg) is minimized at a markedly 
higher u range than for BVS(Zn), making it clear why the ex- 
perimental u(Mg) is conspicuously larger than u(Zn). 

Several authors have used first-principles density functional 
theory (DFT) to predict the relaxed structure of hypothetical 
wurtzite MgO.^^ An earlier DFT study attempted to relax 
the energy-minimized wurtzite MgO structure, but only to ob- 
tain a hexagonal structure corresponding to the limit of c-axis 
compression (c/a = 1.20, u = 0.5).^^ However in later stud- 
ies, the energy minimizations for wurtzite MgO were achieved 
within appropriate boundary conditions; c/a = 1.514, w = 
0.398 by Janotti etal.^^, and c/a = 1.520, u = 0.395 by 
Gopal and SpaldiiPS. Therefore, wurtzite MgO, if it ever oc- 
curs, is expected to have an abnormally small c/a ratio and 
large u value, compared with common wurtzites including 
ZnO. Interestingly, the above computational approximations 
of the hexagonal MgO structure agree well with our PDF anal- 
yses of Zni_;cMg;cO, which underscores that the Mg04 frag- 
ment has smaller c/a and larger u than Zn04. Malashevich 
and Vanderbilt'^'^ and Fan etal^^ performed DFT computa- 
tions of Zni_;cMg;cO supercell models, also deriving consis- 
tent conclusions that c/a decreases with Mg-substitution, and 
that w(Mg) > uiZvi). The results from the neutron PDF refine- 
ment and the above DFT computations are plotted together in 
in which the distinct internal/external geometries of 



Fig. 14 



Mg and Zn were calculated as functions of u (Fig. 13). In 



Mg04 and Zn04 are contrasted. As a general trend, isolated 
Mg04 tetrahedra always have larger u parameters than Zn04 
units. The extremely good agreement between the local coor- 
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■ MgO^ in Zn^MgOg (neutron PDF) 

• ZnO^ in Zn^MgOg (neutron PDF) 

□ IVIgO^ in Zn^^Mg^O supercells (DFT) 
ZnO^ in Zn^^Mg^O supercells (DFT) 

□ MgO^ in hexagonal MgO (DFT) 
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FIG. 14: Aspect ratio cja and internal parameter u for Mg04 (filled 
squares) and Zn04 (filled circles) tetrahedra in Zno.875Mgo.125O, as 
analyzed by neutron PDF refinements. DFT results for wurtzite MgO 
and Zni^^Mgfi supercells are compared together: Open squares cor- 
respond to the Mg04 geometry in hexagonal MgO (i: ref."^^ and ii: 
ref.E2j), and the dotted open squares and dotted open circles represent 
the Mg04 and Zn04 geometries found in Zni_^Mg^O supercells"^^-, 
respectively. The gray line defines the border where the axial cation- 
anion bond along c and three equatorial bonds become equidistant. 



dination of Mg in Zno.875Mgo.125O determined here from PDF 
studies and DFT calculations on end-member wurtzite MgO 
is remarkable. This is in part attributed to our use of the ratio 
cja as a comparison parameter. This ratio has fewer system- 
atic errors than do the values of the lattice parameters taken 
separately. 

The total polarization of polar oxides can be broken 
down into electronic, ionic, and piezoelectric contributions, 
related respectively to the polarizability of component ions, 
ionic arrangements, and lattice parameters. Mg-substitution 
decreases the electronic polarization (Pel) of ZnO, because 
of its greater ionic character, whereas the concomitant c-axis 
compression is expected to increase the piezoelectric polar- 
ization (Ppz). Quantitative estimation of the above two contri- 
butions (Pel, ^pz) requires sophisticated theoretical analyses 
and cannot be provided here for ZnO and Zno.875Mgo.125O. 
Meanwhile, the ionic polarization (Pion) of Zni-^MgxO can 
be approximated by a simple electrostatic model using the 
c-axial positions of Zn, Mg, and O. We calculated of 
Zno.875Mgo.125O and ZnO, assuming point charges {q = +2 
for Zn and Mg, and -2 for O), and using the atomic co- 
ordinates as determined by neutron PDF analysis (see Ap- 
pendix). Separate reference structures were used for ZnO 
and Zno.875Mgo.125O, with the corresponding experimental a 
and c values, but with u = 0.375, so that within isolated 



(Zn,Mg)04 moieties the net centers of the cation and anion 
charges coincide. The calculated P^ values (= Pion here), 
along c, are -5.18jiC/cm^ for ZnO and -5.25jiC/cm^ for 
Zno.875Mgo.125O (APion = -0.07 jiC/cm^). For compari- 
son, Malashevich and Vanderbilt^^ have employed the Berry- 
phase approach, to obtain electric polarization of a Zn5Mg06 
(x = I) supercell. They find APion = -0.22jiC/cm^ 
APei = +0.01 jiC/cm^ and AP^^ = -0.81 jiC/cm^ rela- 
tive to Ps(ZnO) = -3.22 [iC/cm^. In fact, the following re- 
lations hold true in general for other Zni-^^MgjfO supercells 
(x= |, |, ^) studied by them: AP^i > 0, APion < 0, APpz < 0, 
and lAPeil ^ lAPionl < |APpz|. We therefore expect that the 
total polarization of Zni-^^Mgj^O should be larger than that of 
ZnO, with the increases of piezoelectric and ionic polariza- 
tions dominating the decrease of the electronic contribution. 
We note that P^iZnO) from Malashevich and Vande rbilt dif- 
fers somewhat from the results of other DFT studies ^^ l ^^ l ^^^ l 

which reported Ps(ZnO) = -5.0 5.7 jiC/cm^. But since 

we are primarily interested in the relative portions of Pq\, Pion, 
and Ppz in the P^, and their changes with respect to the Mg- 
substitution in ZnO, we ignore the small discrepancies in the 
absolute values. 

The result, an increase of Pion upon Mg-substitution, dif- 
fers from our previous estimates based solely on synchrotron 
x-ray diffraction analyses.^ ^ This inconsistency arises mostly 
from the improved models for the Zni-J^g^O supercell struc- 
tures that are permitted by the superior neutron scattering re- 
sults. Figure[8]clearly shows the greatly enhanced data quality 
of the neutron scattering measurements, compared with x-ray 
experiments. The neutron data lead to better precision (see 
Table|l|andref.^), as well as better accuracy, in the subsequent 
structure analyses. 

It is also noteworthy that the structure modeling for the PDF 
or Rietveld analyses can confine the range of structural in- 
formation that is sought. For example, Rietveld refinement 
of Zni_;cMg;cO used only one positional parameter, and there- 
fore Mg04 and Zn04 tetrahedra are characterized by the same 
c/a and u. In the x-ray PDF refinement, two position param- 
eters were refined, so that we could differentiate u(Mg) and 
w(Zn), but not their c/a ratios. However, for the neutron PDF 
refinement, we used four positional parameters as variables, 
and were able to distinguish the Mg04 and Zn04 geometries 
through both c/a and u. Again, the choice of a structural 
model relies on the quality of the diffraction data, for this ap- 
plication on a system with relatively low-atomic-number ele- 
ments. The neutron data showed sufficient integrity for testing 
various structure models, which were not possible for the x- 
ray data alone. 



IV. CONCLUSION 

Detailed structural analyses of Zni-x^g^O (0 < x < 
0.15) solid solutions, focusing on the dissimilar local geome- 
tries of Mg and Zn, have been performed by using Raman 
scattering, ^^Zn/^^Mg NMR spectroscopy, and neutron Ri- 
etveld/PDF analyses. Line broadenings of the Raman and 
^^Zn/^^Mg NMR spectra imply that Mg-substitution into a 
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ZnO-rich lattice gradually increases the configurational dis- 
order and crystal defects in the lattice. Rietveld refinement 
of time-of-flight neutron diff'raction patterns and Raman peak 
profile analyses show that the macroscopic lattice distortions 
in Zni-x^gxO wurtzites develop in a way such that the c/a 
ratio decreases with Mg content x. A real- space neutron 
PDF analysis using a supercell model of ZnyMgOg reveals 
that the Mg and Zn atoms in Zni-j^Mg^^O have markedly dis- 
tinct local geometries: Mg04 tetrahedra are smaller in height 
and volume and have a larger internal parameter u. The 
wurtzite structural parameters c/a and u for Mg04 deviate 
from their ideal values, in agreement with computational pre- 
dictions of wurtzite MgO structures. Previous DFT studies on 
Zn\-xMgxO supercells and the present neutron PDF analysis 
strongly suggest that Mg- substitution will increase the spon- 
taneous polarization of ZnO. 
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APPENDIX 



Complete list of atomic coordinates for the 16-atom 
wurtzite supercell ZnyMgOg (a = b = 6.5054 A, c = 
5.2015 A, a = |3 = 90°, y = 120°), as determined by the 
neutron PDF analysis. 



Atom 


X 


y 


z 


Type 


Zn 


0.1667 


0.3333 





Znl 




0.6667 


0.3333 





Zn3 




0.1667 


0.8333 





Znl 




0.6667 


0.8333 





Znl 




0.3333 


0.1667 


0.5 


Zn2 




0.8333 


0.1667 


0.5 


Zn2 




0.8333 


0.6667 


0.5 


Zn2 


Mg 


0.3333 


0.6667 


0.486(8) 


Mg 


O 


0.1667 


0.3333 


0.386(3) 


02 




0.6667 


0.3333 


0.3825(8) 


03 




0.1667 


0.8333 


0.386(3) 


02 




0.6667 


0.8333 


0.386(3) 


02 




0.3333 


0.1667 


0.8825(8) 


03 




0.8333 


0.1667 


0.8825(8) 


03 




0.3333 


0.6667 


0.858(5) 


Ol 




0.8333 


0.6667 


0.8825(8) 


03 



